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ABSTRACT 

We present spectra of the afterglow of GRB 021004 taken with WHT ISIS and VLT 
FORSl at three epochs spanning 0.49-6.62 days after the burst. We observe strong 
absorption hkely coming from the host galaxy, alongside absorption in HI, Si IV and 
CIV with blueshifts of up to 2900 km from the explosion centre which we assume 
originates close to the progenitor. We find no significant variability of these spectral 
features. We investigate the origin of the outflowing material and evaluate various 
possible progenitor models. The most plausible explanation is that these result in 
the fossil stellar wind of a highly evolved Wolf-Rayet star. However, ionization from 
the burst itself prevents the existence of HI, Si IV and CIV close to the afterglow 
surface where the fast stellar wind should dominate, and large amounts of blueshifted 
hydrogen are not expected in a Wolf-Rayet star wind. We propose that the Wolf- 
Rayet star wind is enriched by a hydrogen-rich companion, and that the GRB has a 
structured jet geometry in which the gamma rays emerge in a small opening angle 
within the wider opening angle of the cone of the afterglow. This scenario is able to 
explain both the spectral line features and the irregular light curve of this afterglow. 
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1 INTRODUCTION 

The location of Gamma-ray bursts (GRBs) at cosmo- 
logical distances was settled by the measurement of ab- 
sorption lines in optical afterglow spectra, first done for 
GRB 970508 llMetzger et al.lll997l) . Many bursts are now 
followed up spectroscopically, not only for redshift deter- 
mination but also as a probe of conditions in the host 
galaxy, the circumburst environment and intervening ma- 
terial. Spectral features imprinted on the afterglow spec- 
trum can be a useful tool for understanding Gamma-ray 
Burst origins. In several, and at least 2, cases afterglow 
spectroscopy, sometimes in combination with photometry, 
has revealed the presence of a supernova connected with 
the GRB (e.g. GRB 980425 and SN1998bw, Galama et al. 
1998; GRB 011121 and SN2001ke, Garnavich et al. 2003; 
GRB 030329 and SN2003dh, Hjorth et al. 2003; Stanek et al. 
2003 and GRB 031203 and SN20031w, Malesani et al. 2004). 

GRB 021004 was det ected on 2002 October 4, 12:06 UT 
by the HETE-2 sateUite JShirasaki et al.ll2002^ ■ An optical 
afterglow was discovered ten minutes after detection of the 
prompt gamma-ray emission with an R band magnitude of 



15.34 JFox et al.ll2002ah . Intervening absorption systems at 
z = 1.38 and z = 1. 60 were found in the subsequent after- 
glow spectra jFox et al. 200j^, and later spectra showed 
further absorption featu res at redshift z = 2.323 taken to be 
the host galaxy redshift JChornock fc Filippenk J2002^ . The 
afterglow light curve is not simply a power law in shape, but 
shows several bumps, particularly in the well sampled region 
around 1 day after the burst (Fig. 0. The temporal decay 
at early times (before ~ 1 day) follows what is expected 
for a jet exp anding into a stellar w ind medium with a 1/r^ 
dependence jLi fc Chevaliedliool . but when including the 
later times several different explanations are possible. 

GRB 021004 was localised in the optical band very 
early, allowing unusually prompt spectroscopic follow-up. 
The optical afterglow is bright, and rich in line features in- 
cluding Si IV and CIV. It is the first afterglow for which 
the spectrum is dominated by absorption lines spanning a 
3000 km s"^ velocity range, first reported by Salamanca 
et al. (2002). A number of high redshift afterglow spectra 
show absorption features from Si IV and C IV i n their spectra, 
including GRB 000926 dCastro et al.l |2003^ . GRB 011211 
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Figure 1. B and _R band photometry of GRB 021004 compiled 
from tlie literature (Bersier et al. 2003, Fox et al. 2003, Holland 
et al. 2003, Pandey et al. 2003 and Uemura et al. 2003). The cov- 
erage of the WHT and VLT spectra presented here is shown by 
vertical dotted lines. 



l|Holland et al.ll2n02h and GRB 030323 (including an out- 
flowing system with v = 130 ± 60 km s~^, Vreeswijk et al. 
2004). Such lines are also common in the interstellar medium 
(ISM) of many galaxies and al so in a compos ite spectrum 
of z ~ 3 Lyman Break galaxies dShaplev et al.l Eo03'l. but in 
these cases never reach the high velocities observed in the 
afterglow of GRB 021004. Despite many detailed studies of 
afterglow spectra of GRB 021004, the origin of the spectral 
features remains unclear. 

We present early afterglow spectra taken with the 4.2- 
m William Herschel Telescope (WHT) and late-time spectra 
taken with the 8.2-m Very Large Telescope (VLT). Spectra 
were obtained 11.78 hours after the burst - the earliest op- 
tical spectra of GRB 021004 taken by a 4-m class telescope, 
and also at 1.55 and 6.62 days after the burst. We analyse in 
detail the multiple absorption complexes at all three epochs, 
and investigate their origins. 



reduction was performed using IRAF, except for cosmic ray 
removal, where we used the L. A. Cosmic routine to remove 
both point source and irregular-shaped cosmic ray traces 
(van Dokkum 2001). For each night, we combine the indi- 
vidual exposures of GRB 021004 to achieve a higher signal 
to noise. 

The useful remaining wavelength ranges are 3650- 
6000 A for the blue grating and 6500-8500 A for the red 
grating. The wavelength resolution is approximately 2.2 A 
at epoch 1 and 3.5 A at epoch 2. The average signal to noise 
is 12.7 at epoch 1, 7.6 at epoch 2. 

A third epoch observation was made with the Very 
Large Telescope (VLT-Antu) in Paranal on 2002 October 
11, 02:57:09 UT. The observations were made with the Focal 
Reducer/Low Dispersion Spectrograph 1 (FORSl) instru- 
ment in the Longslit Spectroscopy (LSS) mode. Six 1200 s 
exposures were taken through a 0.7 x 416.8 arcsec slit, with 
the 600B grism (3450-5900A). Data reduction was carried 
out using a combination of routines within the IRAF and 
Starlink FIGARO packages. The 6 object images were com- 
bined into a single spectrum. Wavelength calibration was 
perfomed using He and HgCd lamps, and the standard star 
LTT 7987 was used for flux calibration with airmasses rang- 
ing from 1.38-1.57. The resulting spectrum has a signal to 
noise ratio of 4.6 at the longest wavelengths down to ~ 2.5 
at H Lya, and a wavelength resolution of ~ 3.8 A. 

Spectral analysis was done using the Starlink DIPSO 
spectral fltting package. To obtain the correct continuum 
fluxes, the data are compared with B and R band photom- 
etry collected from the literature (Fig. 0. The photomet- 
ric points closest to the observation dates were interpolated 
to give the magnitudes Bcpochi ~ 19.50, iJepoch2 ~ 20.83 
and i?cpochi ~ 18.55 expected in the WHT spectra and 
B ~ 22.44 at the time of the spectroscopic VLT observation. 
The observed fluxes in the WHT spectra are as expected at 
epoch 2, but the earlier spectra are fainter by a factor 2.2 
for the blue grating and 1.7 for the red grating. The flux 
in the VLT spectrum is 1.6 times lower than expected from 
photometric measurements. 



2 OBSERVATIONS 

The WHT spectra were taken on 2002 October 4, 23:52:60 
UT and October 6, 01:23:55 UT with the ISIS spectrograph. 
Overlapping blue (B300B) and red (R316R-HGG495 fllter) 
gratings were used with the EEV12 and Marconi2 CCDs 
respectively. Two exposures totalling 4000 s at epoch 1 and 
4800 s at epoch 2 were taken with each grating, and with 
slits of width 0.87 arcsec at epoch 1 and 1.03 arcsec at epoch 
2. 

Flat flelds were taken during the first night, and 
CuNe-|-CuAr arc lamps were obtained before and after each 
set of target exposures. The standard stars SP0305-I-261 and 
SP2323-I-157 were observed directly preceding or following 
each set of target and arc exposures on both nights. The air- 
mass ranged between 1.01 and 1.09 and conditions were not 
photometric. We applied a correction for atmospheric ex- 
tinction using the mean extinction curve for La Palma, and 
a Galactic extinction correction using E[B ~V) = 0.060 as 
determined by Schlegel, Finkbeiner & Davis (1998). Data 



2.1 Spectral line fitting 

We show the first epoch WHT ISIS spectrum and its error 
spectrum in Fig.|5|with all identified lines labelled. We fit a 
local continuum around each line complex with a first order 
polynomial. Each doublet was modelled with two Gaussian 
profiles, with the line full widths at half maxima (FWHM) of 
the two components tied. For each species, all doublet fits (fit 
parameters of FWHM, central wavelength and peak inten- 
sity) were optimised simultaneously using the least squares 
method. The H Lya region has to be treated differently as 
individual lines are much less clear, so here we fitted a num- 
ber of broad absorption troughs representing some unknown 
number of blended lines, plus the narrow emission line. The 
error on line equivalent widths was measured within the 
DIPSO package using the statistical uncertainties in the con- 
tinuum fit and on the flux values within the line, and also 
the systematic error in contin uum placement and zero l evel 
in the case of the WHT data jHowarth fc PhiUipj[l983) . 
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Table 1. Model fits to the identified spectral line complexes at both epochs. Velocities and equivalent widths are in the restframe of 
the source. EW errors are quoted at the 2cr level and EWs marked with a * indicate lower limits for saturated lines. indicates the 
two heavily blended lines of CIV for which only the total EW should be considered. Velocities and redshifts are given to la accuracy 
and precision of observed wavelengths is limited by the wavelength resolutions specified in Section 2. Epochs 1 and 2 refer to WHT ISIS 
data, epoch 3 refers to VLT FORSl data. 
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Figure 2. The WHT blue grating spectrum at epoch 1 (soUd Hue) and the error spectrum (dashed Une), fluxes corrected according to 
photometric Hght curve in Fig. 1. Identified lines are labelled above. 



3 RESULTS 

The results of spectral fitting are detailed in Table Q We 
observe H Lya, Si IV and CIV absorption features at 3 red- 
shifts, and H Lya and C IV at a further redshift (this is the 
weakest velocity system and is not detected in Si IV or Lya 
at epoch 3), all around z ~ 2.3. Other than the intervening 
systems, all absorption lines have a component at z — 2.327, 
the highest redshift absorption system we observe, which we 
take to be the redshift of the host galaxy. If the other 3 sys- 
tems are a type of outflow (we do not rule out the possibility 
that the velocity system is associated with the GRB pro- 
genitor rather than the host galaxy, but see Section 4.2), 
their observed velocities would be approximately 570, 2400 
and 2800 km s^^ (as reported in Salamanca et al. 2002). 

The lines at we observe at 1190-1197A could be either 
H Lya at 5500 km s"'^, intervening Lya forest lines or Sill 
at km s-i {z = 2.327 ± 0.003). In the former case, we 
might expect to see CIV or Si IV at this velocity as well, 
which we do not; in the latter case, the stronger Sill A1260 
line should also be seen, but it is not. We conclude that the 
nature of the lines at 1190-1197A is uncertain, though most 
likely Lya forest lines, so we claim no H at velocities above 
3000 km s-^ 

3.1 Variability 

We find no significant variability in equivalent width (EW, 
within the 2a errors) between the absorption features com- 



mon to 2 or more epochs, with the exception of AIII A1670 at 
z = 2.327 which we were unable to clearly detect in the 2nd 
and 3rd epochs. Fixing the position of Al II to that found at 
epoch 1 allowed us to fit a line there, and the resulting EW 
is inconsistent with that at epoch 1 within the 2a errors. 
The Al II line is also detected in several other spectra of this 
afterglow (M0ller et al. 2002; Matheson et al. 2003; Mirabal 
et al. 2003; Schaefer et al. 2003; Fiore et al. 2004) taken at 
times both before (in the case of the 2.56-m Nordic Optical 
Telescope (NOT) data) and after our 2nd epoch WHT spec- 
tra, and in each case the measured EW is consistent with 
what we measure in the epoch 1 WHT data. The discrepancy 
may be due to imperfect cosmic ray/bad column removal at 
epoch 2, and/or there may be moderate variability in this 
line. Fell A2344 is another low ionization line which is not 
detected in all optical spectra of GRB 021004, though we 
do find it in the first 2 epochs presented here. The appar- 
ent variability of this line is largely due to the presence of 
a strong skyline near that position. The lack of significant 
variability in the strong lines indicates that the absorbing 
material in our line-of-sight remains in front of the afterglow 
throughout the ~ 6 days spanned by our spectra. 

The EWs we measure are consistent with those of cor- 
responding lines or line complexes reported by M0ller et al. 
(2002), Matheson et al. (2003), Mirabal et al. (2003) and 
Fiore et al. (2004). The resolution of these spectra are such 
that we are not resolving every component of each line blend 
- the line complex at the host galaxy redshift is resolved 
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Figure 3. Outflow velocities of tlie CIV and Si IV absorption 
doublets (velocities plotted are for the longer wavelength compo- 
nents only) and corresponding hydrogen Lya systems, divided by 
the local continua. The solid lines show the WHT observation of 
4/10/02, the dashed lines show the WHT observation of 6/10/02 
and the dotted lines show the VLT observation of 11/10/02. The 
4 velocity systems are denoted by roman numerals, the Lya emis- 
sion line (too large in the VLT data to be shown fully on this plot) 
is labelled 'em'. Intervening C IV is marked with a * and ** marks 
the long-wavelength end of the bluemost absorption complex. 



into 3 components in the VLT UVES spectrum iFiore et al.l 
|20o3)- Unfortunately the cahbration uncertainties in our 
late-time spectrum are too large to detect the colour change 
blueward of Lya of ~ 0.3 magnitudes claimed by Mathe- 
son et al. (2003) and Schaefe r et al. (2003) and supported 
by photometric observations iBersier et al.l 120031) . but not 
found in the HST observations of Fynbo et al. (2004). 

3.2 Hydrogen Lya emission 

We clearly detect the host galaxy ISM through its emis- 
sion in hydrogen Lya. Taking the host galaxy redshift to be 
z = 2.327, the Lya emission line centre then appears red- 
shifted, by an amount that is somewhat different between 
various authors. However, as it lies right on the edge of the 
Lya absorption complex, which probably eats into the emis- 
sion line, its central wavelength will depend on the resolu- 
tion of the spectrum, so this effect is not surprising. Within 
the resolution of these spectra (a few hundred km s~^), its 
velocity is in fact consistent with zero. 

The Lya emission line strength is difficult to estimate 
due to the line blending. However, by fitting in isolation the 
part of the line observed above the continuum level, we can 



measure a lower limit to the flux: 

Fepoch 1 ^ 2.35 ± 0.87 X 10"'^ erg cm'^ s'S Fepoch 2 ^ 
1.66±0.34x 10^^^ erg cm^'^ s~^. These lower limits are con- 
sistent within the quoted Icr errors, suggesting that both the 
emission line and the blended absorption line component at 
the host redshift show no variability, and are co nsistent with 
the value measured in NOT data of this source eretaLl 
l2002h . In epoch 3 we obtain Fopoch 3 1.17x 10"^'' erg cm"^ 
s~^, with a large error of ~ 56 per cent. 

Following the method described in Fynbo et al. (2004, 
and references therein) we calculate a star formation rate 
(SFR) using the Lya emission line flux measured in the 
WHT data. We find a lower limit for the SFR of 6 Mq yr"\ 
with average values of 10 and 7 M© yr~^ at epochs 1 and 
2 respectively, consistent with estimates based on Keck and 
HST data of this afterglow (Djorgovski et al. 2002; Fynbo 
et al. 2004). Our estimate assumes a cosmology with Ho = 
70 km s"^ Mpc"S Q,m = 0.3 and = 0.7. 

From these data we can obtain only a lower limit for 
the column density of HI of TVhi 1? lO^** cm~^, by comparing 
the observed EWs of each line to their theoretical curves of 
growth. This poor constraint is because the HI lines are close 
to or at saturation, they are all heavily blended and cannot 
be well modelled. The column of neutral hydrogen is in fact 
much larger: 1 x 10^** < A^hi (cm^^) < 1.1 x 10^" (M0ller 
et al. 2002; Fynbo et al. 2004). From the unsaturated lines 
of Si IV and CIV (those at the 2 lowest redshifts and part 
of the outflow) we measure columns of log A'^ = 14.0 ± 0.9 
and 14.6 ± 0.9 respectively, where A'^ is given in cm~^, in 
line with previous estimates (Schaefer et al. 2003 C IV only; 
Fiore et al. 2004). 



4 DISCUSSION 

The spectra we have taken provide information on mate- 
rial whose distance from the explosion centre covers a wide 
range. No material can be closer than the emitting surface 
(the blast wave) at the time of the first spectrum, which 
is about 10^^ cm. Some lines are from intervening galaxies 
very far from the host at z = 1.38 and 1.60, and yet oth- 
ers are from the host itself ai z = 2.327. We interpret the 
low-ionization absorption lines and the narrow Lya emis- 
sion line around z ~ 2.33 as coming from within the host, 
but many parsecs from the explosion so that the absorbing 
material was unaffected by the explosion or by the GRB pro- 
genitor. The blueshifted lines we interpret as coming from 
the region affected by the explosion or the progenitor, which 
we shall term the 'high-velocity zone'. The Lya absorption 
could in principle come from either the high- velocity zone or 
the host, but given its velocity structure, identical to that of 
the high- ionization lines of Si IV and CIV, we shall interpret 
it as coming from the high-velocity zone. Ly/3 absorption 
has also been observed with components spanning approxi- 
mately the same velocity range (Mirabal et al 2003.) . Note 

that this implies a very low column of H I in the host in front 
of the immediate star forming region of the GRB progenitor, 
unlike other GRB hosts observed to-date. We discuss these 
outflowing systems and their possible origins in turn below. 
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4.1 Intervening systems, overdense regions and 
partial covering 

There is a possibility tliat tiie 'outflow' may just be several 
systems in the line of sight to the burst. This, however, seems 
highly unlikely to occur so close to the burst redshift (see 
Schaefer et al. 2003; Fiore et al. 2004) and in addition to the 
z = 1.38 and 1.60 systems. Speciflcally, a number of (proto- 
) galaxies within 3000 km s~^ of each other would constitute 
a nonlinearly collapsed overdensity that is very hard to get 
as early as z ~ 2.3 in currently favoured cosmologies either 
in the Hubble flow or within a galaxy cluster including the 
GRB host. 

Let us now consider a clumpy stellar wind environment 
for GRB Q21QQ4. Clumps of material may comprise the stel- 
lar wind - Cherepashchuk (1990) suggested that up to 80 per 
cent of a Wolf-Rayet stars' mass can be lost in clumps in the 
wind. In principle the lines we see could either be due to a 
wind with optical depth less than 1 absorbing a fraction of 
the light more or less uniformly throughout the wind mate- 
rial, or by a partial covering of the emitting surface by opti- 
cally thick clouds. In the latter case, the spectrum should be 
the sum of a featureless continuum from regions uncovered 
by cloud, and a spectrum with saturated, flat-bottomed lines 
from the regions covered by clouds that are optically thick 
in the lines but thin in the continuum. This implies that the 
summed spectrum should have flat-bottomed lines whose 
flux level in the centre is non-zero, when viewed with high 
enough resolution spectroscopy. Whilst our data do not al- 
low us to address this distinction, the UVES high resolution 
afterglow spectra of GRB 021004 (Fiore et al. 2004; Castro- 
Tirado et al., in preparation) show flat-bottomed, saturated 
lines with no flux in the centre, effectively ruling out a par- 
tial covering model. 

We now consider the exotic possibility that the high- 
velocity absorption features are from a number of supernova 
shells produced by an earlier stellar death in the same cluster 
that produced the progenitor of the GRB. It is plausible that 
an unusually rich cluster (e.g. the central cluster of the S Dor 
complex in the Large Magellanic Cloud, LMC) would have 
produced a number of stellar explosions in the past Myr, 
of which the supernova remnants are now expanded beyond 
the cluster core size, so that our line of sight to the next ex- 
plosion will only pass through the approaching sides of the 
shells and thus give rise to blueshifted features. The hydro- 
gen richness is explained in this case, since the ejecta could 
contain hydrogen, or sweep it up. The problem is that the 
velocities are low enough that supernova shells must have 
significantly decelerated. In so doing, they will shock and be 
heated to X-ray temperatures, which makes species like H I, 
Si IV and CIV rather unlikely to be present in those shells. 
The shell geometry has also been seen around Wolf-Rayet 
stars as shell nebulae. Mirabal et al. (2003) demonstrated 
that optical absorption lines like those seen in GRB 021004 
can be produced by a fragmented shell nebula, under very 
particular conditions. The shell nebula would need to have 
been located at a distance from the GRB progenitor such 
that it received sufficient radiative acceleration from the 
burst to achieve the high velocities, but was not fully ionized 
by the GRB flux. All things considered, we do not favour 
these options. 



4.2 A single star as progenitor 

The most natural explanation for the blueshifted absorption 
features is that they arise in a stellar wind. A l/r^ density 
profile is inferred from the light curve at very early times 
(Li & Chevalier 2003) suggesting a stellar wind circumburst 
medium. A GRB progenitor must be a massive, evolved star. 
The highest velocity system observed is an outflow of ~ 
2800 km s~^, which could be the terminal velocity of a stel- 
lar wind and in this case implies the star is a Wolf-Rayet 
(WR) star (Mirabal et al. 2003; Schaefer et al. 2003). Such 
stars evolve to type lb supernovae and are characterised by 
large mass loss rates in high velocity, optically thick stel- 
lar winds. The 3 outflow velocities observed in this after- 
glow could correspond to the fast WR wind (2800 km s~^), 
a slower wind from an earlier phase in the stellar lifetime 
(2400 km s"^) and the slow bubble moving into the ISM 
or mixing of winds from different phases (570 km s~^) as 
shown in simulations of massive star evolution by Van Marie, 
Langer & Garcia-Segura (2005). Massive star winds are pre- 
dicted to show a higher optical depth in th e CIV line than 
in the Si IV line (iLeitherer fc Lamerslll99'lh which seems to 
be true for these data. Line-locking is a characteristic of 
WR wind spectra, since the winds are radiatively driven 
(e.g. Puis et al . 1996), and this ap pears to be present in the 
WHT spectra \M0\\ er et al ]|2002l) . seen as an overlap of the 
CIV doublet lines at the two highest redshifts. If, however, 
the only signature of the radiatively driven WR wind is the 
fastest absorption component, then line-locking would not 
be expected, and would be a coincidence in this case. 

There are two basic challenges to interpreting these lines 
as originating in the stellar wind. The first is that the species 
causing the lines must survive the initial blast of UV to 
gamma rays without becoming fully ionized. The distance 
out to which the material might be fully ionized can be es- 
timated in two different ways, depending on the recombina- 
tion time-scale. If the recombination time is short compared 
to th e 100 s duration of the ionizing flash ("Shirasa ki et al.l 
l2002h . we can use ionization equilibrium as characterised by 
the ionization parameter ^ = L/nr'^ . 

If on the other hand the recombination time is long, 
then a given atom with cross-section a for ionization will 
simply absorb at least one ionizing photon: 

^ 1 (1) 



0"i 



47rr2 < Ei> 



This is the appropriate limit to use since a rough calculation 
of the recombination time for hydrogen gives a time-scale of 
the order of 10^ years. This gives minimal distances for the 
survival of Si IV, CIV and HI of 5x10^**, 2x10^^ and 9x10^" 
cm respectively (using photoionization cross-sections and 
threshold energies from Verner et al. 1996 and E^ot ~ £iso)- 
If we assume that these 3 species are co-located (because 
they have the same velocity structure) then the minimum 
distance out to which we could find them if they were ion- 
ized by the burst flux is far beyond the few pc radius of a 
typical stellar wind bubble. This excludes a stellar wind ori- 
gin for these absorption lines if they are intercepted by the 
GRB flux. 

The second challenge is the presence of hydrogen in the 
spectra with approximately the same velocity structure as 
Si IV and CIV, strongly suggesting that significant amounts 
of hydrogen are present in the stellar wind, with a column 
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of neutral hydrogen ^10^* cm^^. This is not usual for a WR 
wind. A WR star will once have been hydrogen-rich, and so 
the outer wind can have hydrogen, but its velocity structure 
will not match that of the later, hydrogen-deficient wind. 

WR stars in our own galaxy are hydrogen deficient 
and in most cases H is not detected at all. However in the 
SMC a large fraction, if not all single WN-type WR stars 
nitrogen-rich WR sta rs) show hydrogen in their spectra 
Foellmi. Moffat fc Guer rero 2 fl03aV This was unexpected 
given the lack of H in Galactic WN-types, leading to the 
suggestion that the WN population in the SMC is funda- 
mentally different from that of our own galaxy. So in the 
SMC, where the metallicity is low causing less severe mass 
loss and fewer WC-type (carbon-rich) WR stars, massive 
stars can remain H-rich in the WR phase. GRB host galaxies 
are thought to be metal-poor like the SMC, but the strength 
of C IV and Si IV and non-detection of nitrogen (particularly 
NIV, NV) in the afterglow spectra of GRB 021004 is incon- 
sistent with a WN-type WR star. In fact carbon produces 
the strongest absorption features, suggesting we are observ- 
ing the wind of a WC-type WR star (e.g. Mirabal et al. 
2003; Schaefer et al. 2003), thought to be the last evolu- 
tionary stage of the WR phase. On the other hand we do 
not detect CIII A1247 which is strongly detected in samples 
of WC-types with e.g. F USE dWillis et al.ll2004) and lUE 
iNussbaumer et al .''l982Y 

There is no significant variability in the absorption line 
EWs between our very early spectrum at 0.49 days after 
burst, and those beginning 1.55 and 6.62 days after burst 
(besides the AIII line discussed in Section 3.2). The three 
epochs appear to be in different 'phases' of the peculiar af- 
terglow light curve (Fig.0. The jumps in the light curve at 
0.8-1 and 2.5-3 days cannot then be shells through which 
the afterglow light passes, because the same absorption is 
seen at all epochs. Unless the shells are located very far 
from the star so that the afterglow never overtakes them, 
which then makes it difficult to produce the fast veloci- 
ties we observe. 21-cm measure ments have shown HI shells 
around WR stars, e.g. WR 102 iGosachinskii fc Lozinskaval 
|2002|) . but these have expansion velocities of ~ 50 km s~^. 
The Galactic WR star WR 3 may also have an unusually 
high H content if it is a single star rather than in binary 
with an O star, but again o nly rather low velocitie s of up 
to 400 km s~ are observed jMarchenko et al.ll200l) . never 
reaching the 2400-2900 km s~^ velocities we observe. 

4.3 A binary system progenitor 

So, whilst a WR star is the obvious choice to explain the 
blue shifted, line-locked absorption complexes, this picture is 
difficult to reconcile with the large amount of neutral hydro- 
gen observed, particularly at very high velocities. This co- 
nundrum leaves us with few options. Either distant carbon- 
rich WR stars are very different to those in our own galaxy 
and are able to retain a large fraction of their hydrogen, or 
the hydrogen originates elsewhere - perhaps in a companion 
star wind. There is the possibility that we are seeing a WR 
star in a close binary with a hydrogen-rich main-sequence 
star. Even for an O star companion, the momentum loss 
rate Mvoo of the WR star far exceeds that of the compan- 
ion, and so the wind velocity structure at distances much 
greater than the binary separation will be dominated by the 



WR star. However, a few to 10 per cent of the mass will come 
from the companion and so is hydrogen rich. When mixed 
into the WR wind it can cause the observed phenomenon of 
H LyQ emission accompanying otherwise typical WR lines. 
The a priori probability that a WR star is in a binary with 
an O star is still a debatable number. Foellmi et al. (2003a) 
give theoretical binary frequency lower limits for WR-l-0 
systems in the Magellanic Clouds of 0.41 ± 0.13 (LMC) and 
0.98 ± 0.32 (SMC), assuming a mass loss - metallicity rela- 
tion of M oc Z"'^ and excluding rotation effects. They find, 
however, significantly lower observed fractions of 30 and 40 
per cent respectively (Foellmi et al. 2003a, b). Various mea- 
sured values point towards the conclusion that the binary 
frequency of WR stars is identical to that of their progeni- 
tors and independent of the metallicity. 

There is also the possibility that the WR star is in a 
binary with a less massive star that is providing the observed 
hydrogen. The probability of a WR star having a lower- 
mass companion (perhaps 1-3 M©) is even less well known, 
since such low-mass, low-luminosity companions are hard to 
detect. The only evidence for the existence of such binaries 
is the population of soft X-ray transients in our Galaxy, 
which consist of a black hole orbited by a low-mass star 
in a very tight orbit. Their ancestors must have been WR 
stars with a low-mass main-sequence binary companion. It 
has been suggested that these s ystems produce a GRB when 
the black hole forms iLee. Brown fc W iicrs 2002). 

One necessary consequence of this H-admixing model is 
that probably only within some angle from the orbital plane 
of the binary will H be mixed into the outflowing wind. 
In the polar direction pure WR composition will dominate. 
This will constrain our viewing angle to the system to being 
not too far out of the orbital plane. For an O star, still 
having a relatively strong wind, this is not severe, but for 
a low-mass star with a very weak wind, the system would 
have to be viewed at or close to edge-on. This is unlikely if, 
as predicted by the coUapsar model, the jets emerge from 
the polar axis of the progenitor star. 

4.4 A structured jet model 

The ionization problem remains however, and it seems that 
in order to observe fast-moving species such as Si IV and C IV 
and even HI, one must be observing a region undisturbed 
by the gamma ray jet. In order to achieve this we suggest a 
structured jet geometry, with one narrow jet (< 3°) and one 
slower, much wider jet of material which are both visible to 
us (Fig. 0J. In fact, the afterglow has a surface area, and 
so there are several lines-of-sight along which the light will 
reach us. Gamma rays from the inner jet will ionize material 
in its path. But if the inner jet opening angle is small then 
there will also be a line-of-sight along which we see material 
unaffected by high energy radiation. We propose that there 
is a second wider cone of slower jet material surrounding the 
narrow gamma ray jet. This wide jet is also seen through the 
absorber, probably a fast stellar wind within 2 older more 
slowly moving wind phases, which produces the absorption 
lines we observe. 

The power law decay of afterglow light curves can be 
very smooth, as in the cases of GRB 020813 (e.g. Laursen & 
Stanek 2003) and GRB 990510 (e.g. Stanek et al. 1999). But 
some show considerable deviations from this with many ir- 
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Figure 4. A sketch of the structured jet geometry within the 
stellar wind, where our Unes-of-sight (dashed lines) intersect both 
fully ionized material infront of the narrow, fast jet (dark grey) 
and material ionized to a much lesser degree by the the slower jet 
(light grey). 



regularities, as seen in GRB 011211 (e.g. Holland et al. 2002) 
and GRB 030329 (e.g. Lipkin et al. 2004) and this is the case 
for this burst (Fig^ . The structured jet scenario we describe 
is consistent with the observed light curve of the afterglow 
of GRB 021004 (Van der Horst, private communication). We 
also note that the general decay of the light curve may be 
adequately fit with other models including several energ y 
injections llBiornss oiij_Gudrnu^^ . 
patchy shell models jNakar^ Pirai^^^^no?200?rand den- 
sity e nhancements in the ISM ( Lazzati et al. 2002) or in a 
wind teevl fc Pernall2003D . 

From the light curve alone we cannot pinpoint the cause 
of the irregularities in the temporal decay. Considering the 
light curve in combination with the spectra provides further 
constraints and leads us to prefer the structured jet idea, 
which can explain both types of observation and does not 
require the absorbers to be at a large distance from the pro- 
genitor to produce the absorption lines seen in our spectra. 



5 CONCLUSIONS 

We have observed the optical afterglow of GRB 021004 spec- 
troscopically at three epochs ranging from 0.49 to 6.62 days 
after the GRB, when the emitting surface is of order 10^^ cm 
away from the explosion centre. The most unique feature in 
the spectra is the fact that lines of H I, Si IV and C IV are ob- 
served at blueshifts of up to 2900 km s~^. While no model 
explored by ourselves or earlier authors explains this phe- 
nomenon in a very simple manner, we find that by far the 
most plausible explanation is that these lines result in the 
fossil stellar wind of a highly evolved Wolf-Rayet star. It 



is unusual that hydrogen would be present in such a wind. 
This could be explained if WR stars do remain hydrogen- 
rich up to explosion at high redshift and/or low metallicity, 
or if a nearby companion star had enriched the WR wind 
with hydrogen. In any case, our data seem to further en- 
hance the case for the origin of GRBs in the core collapse 
of very massive stars. Ionization from the burst itself pre- 
vents the existence of H I, Si IV and C IV close to the after- 
glow surface where the fast stellar wind should dominate. 
A structured jet geometry is a promising explanation: the 
gamma rays emerge in a small opening angle leaving a wider 
cone of slower, less energetic material also along our line-of- 
sight. This wider jet does not fully ionize the material in its 
path, and can give rise to the absorption lines we observe. 
This model is able to explain both the spectral features and 
reproduce the light curve shape. 
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